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Abstract 

A system of N unidimensional global coupled maps (GCM), which support multiattractors is 
studied. We analize the phase diagram and some special features of the transitions (volumen ratios 
and characteristic exponents), by controlling the number of elements of the initial partition that 
are in each basin of attraction. It was found important differences with widely known coupled 
systems with a single attractor. 
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I. INTRODUCTION 



The emergence of non trivial collective behaviour in multidimensional systems has been 
analized in the last years by many authors [[! ] |2j . Those important class of systems are 
the ones that present global interactions. 

A basic model extensively analized by Kaneko is an unidimensional array of N elements: 

e N 

X t+1 (k) = (1 - c) f(X t (k)) + f{X t (l)) (1) 

iV i=i 

where k = 1, . . . , N, is an index identifying the elements of the array, t a temporal discret 
variable, e is the coupling parameter and f(x) describes the local dynamic and taken as the 
logistic map. In this work, we consider f(x) as a cubic map given by: 

f(x) = (l-a)x + ax 3 (2) 

where a G [0,4] is a control parameter and x G [—1,1]. The map dynamic has been 
extensively studied by Testa et.al. @], and many applications come up from artificial neural 
networks where the cubic map, as local dynamic, is taken into account for modelizing an 
associative memory system. Ishi et. al. || proposed a GCM model to modelize this system 
optimazing the Hopfield's model. 

The subarmonic cascade, showed on fig-| prove the coexistence of two equal volume stable 
attractors. The later is verified even as the GCM given by Eq.^Jhas e > 0. Janosi et. al. || 
studied a globally coupled multiattractor quartic map with different volume basin attractors, 
which is as simple second iterate of the map proposed by Kaneko, emphazasing their analysis 
on the control parameter of the local dynamic. They showed that for these systems the mean 
field dynamic is controlled by the number of elements in the initial partition of each basin 
of attraction. This behaviour is also present in the map used in this work. 

In order to study the coherent-ordered phase transition of the Kaneko's GCM model, 
Cerdeira et. al. [fjj analized the mechanism of the on-off intermitency appearing in the 
onset of this transition. Since the cubic map is characterized by a dynamic with multiple 
attractors, the first step to determine the differences with the well known cuadratic map 
given by Kaneko is to obtain the phase diagram of Eq.|l| and to study the the coherent- 
ordered dynamical transition for a fixed value of the control parameter a. The later is done 
near an internal crisis of the cubic map, as a function of the number of elements Ni with 
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initial conditions in one basin and the values of the coupling parameter e, setting N equal 
to 256. After that, the existence of an inverse period doubling bifurcation as function of e 
and Ni is analized. 
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FIG. 1: The subharmonic cascade, we graphic 256 iterations of the map versus a for each basin, 
after a transient of 5000 steps. 

II. PHASE DIAGRAM 

The dynamical analysis process breaks the phase space in sets formed by synchronized 
elements which are called clusters. This is so, even when, there are identical interactions 
between identical elements. The system is labeled as 1-cluster, 2-cluster, etc. state if the 
Xi values fall into one, two or more sets of synchronized elements of the phase space. Two 
different elements i and j belong to the same cluster within a precision 5 (we consider 
5 = 1CT 6 ) only if 

\Xi(t) - XjW < 6 (3) 
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Thus the system of Eq.|], shows the existence of different phases with clustering (coherent, 
ordered, partially ordered, turbulent). This phenomena appearing in GCM was studied by 
Kaneko for logistic coupled maps when the control and coupling parameters vary. A rough 
phase diagram for an array of 256 elements is determined for the number of clusters calcu- 
lated from 500 randomly sets of initial conditions within the precision specified above. This 
diagram displayed in fig-Q, was obtained following the criteria established by this author. 
Therefore, the K number of clusters and the number of elements that build them are relevant 
magnitudes to characterize the system behaviour. 
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FIG. 2: Phase diagram as function of e and a. Number such (2, 4) represents dominant cluster 
number (that is with basin volume ratio more than 10%) 



III. PHASE TRANSITION 

In order to study phase transition, the two greatest Lyapunov exponents are shown in 
fig44] and fig45|. They are depicted for a=3.34 as a function of e and for three different values 
of initial elements N±. 

In the coherent phase, as soon as e decrease, the maximum Lyapunov exponent changes 
steeply from a positive to a negative value when the two cluster state is reached. A sudden 
change in the attractor phase space occurs for a critical value of the coupling parameter e c 
in the analysis of the transition from two to one cluster state. Besides that, in the same 
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transition for the same e c , a metastable transient state of two cluster to one cluster chaotic 
state is observed, due to the existence of an unstable orbit inside of the chaotic basin of 
attraction, as is shown in fig-|] The characteristic time T t in which the system is entertained 
in the metastable transient is depicted in Fig-|B|, for values of e near and above e c . 
For a given set of initial conditions, it is possible to fit this transient as: 

Tt = (—)-" (4) 

This fitting exponent (3, depends upon the number of elements with initial conditions in 
each basin as is shown in the next table for three N± values and setting N = 256. 



iVi 


e c 


P 


128 


0.471829 


0.792734 


95 


0.3697115 


0.606751 


64 


0.3198161 


0.519833 



It is worth noting from the table that f3 increases with N\ up to N/2, and for N\ due to 
the basins symmetry. 

IV. INVERSE PERIOD DOUBLING BIFURCATION 

In order to analize the existence of period doubling bifurcations, the maxima Lyapunov 
exponent \ m ax is calculated as function of N\ and e. For each N\, critical values of the 
coupling parameter, called e^/, are observed when a negative X max reaches a zero value 
without changing sign. This behaviour is related to inverse period doubling bifurcations of 
the GCM. Fitting all these critical pair of values (e^/, iVi), a rough N\ vs e^f graph is shown 
in fig-^, and different curves appears as boundary regions of the parameter space where the 
system displays 2 n {n = 1, 2, 3) periods states . This is obtained without taking into accout 
the number of final clusters. It is clear that greater values of Ni, correspond to smaller 
€bif for the occurrence of the bifurcation. Evidence of period 16 appears for values of N\ 
smaller than 30. In fig-(7] T=2 (symmetric) means period two orbit, with clusters oscillating 
with equal amplitud around zero, T=2 (asymmetric) means period two orbit, with clusters 
oscillating with different amplitud. 
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FIG. 3: X vs. t for two elements of differents clusters near the metastable transition, for a = 3.34, 
N x = 95 and e = 0.3697126. 

V. CONCLUDING REMARKS 

The study of systems with coexistence of multiple attractors gives a much richer dynamics 
and a new control parameter must necessarily be added. Although the dimensionality in the 
parameter space is increased by one, the dynamics is rather simple to characterize. Some of 
the relevant aspects of this kind of systems are shown in this work. The phase diagram that 
was obtained shows the existence of similar phases to those using the cuadratic and quartic 
map, this behaviour suggests some kind of universality in the dynamics of the GCM. Another 
interesting issue found, concerns the metastable transition between two to one cluster state, 
along with a sudden jump in the maximum Lyapunov exponent, as it was displayed in fig.0. 
The characteristic time given by Eq.[| also correspond to the above transition where the 
critical exponent (3 and the critical coupling parameter e shows a strong dependence on the 
number of initial elements in each basin. An inverse bifurcation cascade appears when the 
system is in two or more clusters state where e and N\ are the critical parameters of the 
bifurcation, which means the maximum Lyapunov exponent is equal to zero. 
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FIG. 4: Maximum Lyapunov exponents for three different values of N\. The solid line corresponds 
to Ni = 128, the dashed line to Ni = 95, and the doted line to Ni = 64. 
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FIG. 5: Second Lyapunov exponents for three different values of N% The solid line corresponds to 
Nx = 128, the dashed line to Ni = 95, and the doted line to Ni = 64. 
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FIG. 6: Characteristic time Tt for N± = 128, doted line, and its fitting with the Eq.||, solid line. 
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FIG. 7: Bifurcation diagram for a = 3.34 up to period 8, obtained from 1000 sets of initial 
conditions. 
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